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Abstract 

This note relates topics in statistical mechanics, graph theory and combinatorics, lattice 
quantum field theory, super quantum mechanics and string theory. We give a precise re- 
lation between the dimer model on a graph embedded on a torus and the massless free 
Majorana fermion living on the same lattice. A loop expansion of the fermion determinant 
is performed, where the loops turn out to be compositions of two perfect matchings. These 
loop states are sorted into co-chain groups using categorification techniques similar to the 
ones used for categorifying knot polynomials. The Euler characteristic of the resulting co- 
chain complex recovers the Newton polynomial of the dimer model. We re-interpret this 
system as supersymmetric quantum mechanics, where configurations with vanishing net 
winding number form the groimd states. Finally, we make use of the quiver gauge theory - 
dimer model correspondence to obtain an interpretation of the loops in terms of the physics 
of D-branes probing a toric Calabi-Yau singularity. 
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1 Introduction 



The dimer model is concerned with the statistical mechanics of close packed dimer arrange- 
ments on a bipartite graph. The physical system which can be thought of as a real world 
representation of the dimer model is the adsorption of diatomic molecules on a crystal sur- 
face. 

In the 1960s, the question of how many perfect matchings exist on a plane graph was 
solved independently by Kastele}m 111121, and Temperley and Fisher [|3H1|: the total number 
is given by the Pfaffian of a signed, weighted adjacency matrix of the graph (the Kasteleyn 
matrix). Much of the original interest in the dimer model arose because it provides a simple 
and elegant solution for the 2-dimensional Ising model |5|. 

The problem of enumerating perfect matchings is of course a classical problem in graph 
theory and combinatorics, see e.g. [6J. It can also be phrased in terms of domino tilings [7J. 
During the last years, the interest in the dimer model was revived thanks its manifold con- 
nections to other branches of mathematics and physics. The dimer model is related to 

• configurations of a melting crystal corner and the topological string A-model HHH, 

• real algebraic geometry IITOl[TT| , 

• BPS black holes from D-branes wrapping collapsed cycles |IT2| . 

Furthermore, a correspondence between the dimer model and quiver gauge theories arising 
from D3-branes probing a singular toric surface was discovered and worked out in great 
detail ||T3l[T4l[T5l[T6l[T7| . An explanation of this correspondence via mirror symmetry was 
given in ||T8| . 

That the dimer model must be related to the free massless fermion on the same lattice 
has been known for a long time already. Its relation to the Ising model, which in turn cor- 
responds to a free fermion model clearly indicates this. Moreover, Kasteleyn also showed 
that the partition function of the dimer model living on a graph embeddable on a genus g 
Riemarm surface is given by a linear combination of 2^^ Pfaffians. This is, of course, reminis- 
cent of the number of spin structures on a Riemarm surface. In fact, a one-to-one correspon- 
dence between equivalence classes of Kastele}m orientations and spin structures was proved 
in in. 

In the following, we will make the connection between the dimer model and the free 
fermion precise by identifying the Dirac operator with the Kasteleyn matrix and relating the 
expansion of the determinant appearing in the free fermion path integral to loop configura- 
tions composed of two perfect matchings. This picture of free fermions as a loop gas already 
appeared in a different context, for example in [20] . 

The perfect matchings form a basis for the space of fermion loop configurations. If the 
loop gas lives on a surface of non-trivial homotopy, non-contractible loops of non-zero 
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winding number can occur, and the Pfaffian of the Kasteleyn matrix becomes a polyno- 
mial. Drawing inspiration from the categorification programme of Khovanov l|2Tl l22l [23l , 
we interpret this polynomial as the Euler characteristic of a co-chain complex. The loop 
configurations of the free fermion model can be naturally classified by the weight of the 
constituent matchings of the loop state and used to generate Abelian groups. For the Euler 
characteristic to coincide with the Newton polynomial, a differential operator must be con- 
structed accordingly. In our fermion loop gas, all loop states are paired since they come in 
two opposite orientations, except for states consisting only of double line perfect matchings. 
This is reminiscent of the supersymmetric ground states in supersymmetric quantum me- 
chanics. Indeed, it is possible to re-interpret the fermion loop gas as SQM and to construct 
a differential operator which can be identified with the Q operator. Like this, the Newton 
polynomial becomes the generalized Witten index of the system. Since the fermion loop 
states are bilinear in the perfect matchings, they can also be interpreted as maps from one 
matching to another. 

Thanks to the correspondence between the dimer model and quiver gauge theory, we 
know that the perfect matchings, i.e. our zero modes, parametrize the toric surface being 
probed by the D-branes. In the following, we make use of this correspondence to give the 
loops states an interpretation in terms of the quiver gauge theory. Loops of vanishing net 
winding can be related to sequences of double Seiberg dualities on the associated acyclic 
quivers. In the special case of the hexagon graph, they correspond to certain BPS black hole 
configurations which are parametrized by melting crystal configurations |[T2|. 

The plan of this paper is as follows. In Section |2| we review the necessary background 
and definitions of the dimer model. In Section |3.1[ we relate the Kasteleyn matrix to the 



discretized Dirac operator of the free fermion. In Section 3.2 we go to the limit of vanishing 
lattice spacing in which we recover the partition function of the Dirac fermion, which in the 
light of the interpretation of the dimer model as a free fermion model is no longer surprising. 



In Section 3.3 we expand the determinant arising in the free fermion path integral into cyclic 



permutations, which correspond to closed loops on the graph. Furthermore, we give an 



interpretation of the loops as fermionic states in Subsection 3.4 

In Section|4| the basic idea of categorification is explained, and in Section 4.1 the fermion 
loop system is re-interpreted as supersymmetric quantum mechanics. The simpler case of 
a square graph on the cylinder is worked out first, together with a simple explicit example 



in Section 4.2 In Section [43} we interpret the loops as operators mapping from one perfect 



matching to another. Finally, Section 4.4 gives the full construction on the torus. 

In Section |5| the basics of quiver gauge theories and the correspondence to the dimer 
model is reviewed. The interpretation of the loop states in terms of the physics of D-branes 



probing a singular toric CY is attempted in Section |6j In Section 6.1 the relation between 



perfect matchings and corresponding acyclic quiver graphs is clarified. In Section 6.2 we 
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show that sequences of double Seiberg dualities correspond to cyclic permutations of the 
nodes of the acyclic quiver and to isomorphisms of the full gauge quiver. With these results, 
we can interpret the loops with vanishing overall winding in Section 6.3 The loops with 



non-trivial winding are discussed in Section 6.4 



The results of this paper are summarized in a dictionary given in Section where we 
close with some concluding remarks. 



2 Dimer Model Preliminaries 

Take a plane bipartite graph i.e. one in which all vertices can be colored black and white, 
such that each black vertex is only connected by links to white vertices and vice versa. Let M 
be a subset of the set E of edges of W. M is called a matching, if its elements are links and no 
two of them are adjacent. If every vertex of is saturated under M, the matching is called 
perfect. Such a link which joins a black and a white vertex is also called a dimer. The dimer 
model describes the statistical mechanics of a system of random perfect matchings. In the 
simplest case, we ask for the number of close packed dimer configurations, i.e. the number 
of perfect matchings. 

Let us label the vertices of the underlying graph W on which the dimer model lives con- 
secutively. The topology of $f can be encoded in its adjacency matrix A, where Ax,y = 1 
if the vertices x and y are joined by an edge, and Ax,y = otherwise. The number of per- 
fect matchings is given by the number of ways in which the vertices can be partitioned into 
adjacent pairs. 

The Hafnian of a symmetric matrix A, introduced by Caianiello f2A], is defined as 

Hf(A) = ^7r(l),7r(2)^7T(3),7r(4) • • • ^TT(n-l),n(n)' (1) 

where S is the set of all permutations tt of { 1, ... n } satisfying 7r(l) < tc{3) < ■ ■ ■ < n{n — 
1) and n{2i — 1) < rc{2i) for 1 < i < n, n even. If we take A to be the adjacency matrix of ^, 
Hf (A) counts the perfect matchings, since terms in ^ which contain non-adjacent pairs are 
zero. The permanent of an n x n matrix A is, like the determinant, given by the sum over all 
permutations, except that the sign of the permutations is not taken into account: 

n 

perm(A) = H'^rMiy (2) 

While the determinant can be evaluated in polynomial time via Gaussian elimination, com- 
puting the permanent is #P-complete 1251 . The Hafnian can be regarded as a generalization 
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of the permanent. In fact ||26t . 



Hf (A) = perm(B) for A of the form A = 

This identity accounts for the fact that sometimes, the number of perfect matchings is given 
as the permanent of the adjacency matrix in the literature. Like the permanent, also the 
Hafnian is difficult to evaluate and does not satisfy any useful identities. 

The counterparts of the Hafnian and permanent are the Pfaffian and determinant, which 
take into account the sign of the permutations: 

Pf(A) = ^ sign(7r) A^(i)^^(2)A^(3)^^(4) . . . A„ (4) 

n 

det(A) = X] sign(7r)nA,7r(0' (5) 

where the permutations contributing to the Pfaffian are subject to the same conditions as for 
the Hafnian. Here, we have 

Pf(A') = (-l)(2)det(B) for A' of the form A' = 

Moreover, for every anti-symmetric matrix of even order, the relation 

detA=(PfA)^ (7) 

holds. Already in 1913, Polya suggested to change some signs in a matrix A, such that the 
determinant of the new matrix A would be the permanent of A [27]. Kastele5m |jTll3 actually 
does exactly this: he introduces an orientation on which leads to a signed adjacency 
matrix K, now called the Kasteleyn matrix. The Pfaffian of K gives the number of perfect 
matchings. 

A Kasteleyn orientation fulfills the following condition: the product of all edge weights 
around a face must equal —1 if the number of edges around the face is mod 4. If the num- 
ber of edges equals 2 mod 4, the product must equal 1 ||6|. One can choose an orientation 
by consistently assigning arrows to the edges of the graph, as originally suggested by Kaste- 
lejm lUlll. For our purposes, it turns out to be more convenient to allow the roots of 1 as 
(complex) edge weights. 

Already Caianiello |24| remarked that the expectation value of any product of free Fermi 
fields is a Pfaffian, while the expectation value of any product of free Bose fields is a Hafnian. 
In this sense Kasteleyn's method of attaching signs to the edges of a graph ^ could be seen 
as a fermionization of a bosonic system. Given the correspondence of the dimer model to 



(3) 
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the free fermion on , one could look for its bosonic equivalent, but this goes beyond the 
scope of the present work. 

In the following, we will consider graphs embedded on a cylinder or a torus. The treat- 
ment can be straightforwardly generalized to any genus g Riemann surface. In essence, all 
of the above remains true. The only change is that on the torus (cylinder), we have two (one) 
non-trivial cycles, which we will denote by z and w. In the case of the plane graph, the edge 
weights originated solely from the Kasteleyn orientation. More generally, the Kuperberg 
flatness condition must be respected Il28l . We choose a positive direction on the dimers, say 
• — > o. Now we assign the weight z (vo) to each edge which crosses the cycle z (w) in posi- 
tive direction and the weight 1/z (1/w) to each edge which crosses it in negative direction. 
While the Pfaffian of the Kasteleyn matrix yielded a number in the case of the plane graph, 
it now becomes a polynomial in z and w, the so-called characteristic polynomial or Newton 
polynomial of the graph. The coefficient of each monomial z^'w'? gives the number of match- 
ings with weight (z, w) = {p, q). These are matchings with the number of dimers crossing z 
in positive direction minus the number of dimers crossing z in negative direction equal to p 
(analogous for q). In the literature, what we call the weight is usually referred to as the slope 
of a height function defined on the composition of two matchings. The height function is 
defined as follows. Choose a reference matching PMq. To find the slope of a matching PM, 
compose it with the reference matching, PM — PMq, where the minus serves to change the 
orientation of PMq to o •. This results in closed loops (composition cycles) and double 
line dimers. The rule is that when an edge in PM belonging to a closed loop is crossed such 
that the black node is to its left (right), the height changes by +1 ( — 1). If an edge belonging 
to PMq is crossed, the signs are reversed. This height function is defined up to the choice of 
the reference matching PMq. Crossing the boundary of the fundamental region of the torus, 
this function can jump. If the height function jumps by p units crossing z, it is associated to 
the power zP in the Newton polynomial of the graph (and equivalently for w). Choosing a 
different reference matching results in a common prefactor of zP^w'"' for all monomials. Our 
method of assigning weights to a matching corresponds to choosing a reference matching of 
weight (0,0) that does not intersect the z or w cycle. The matching shown in Figure [T] has 
weight (1,0), where 1 = 2 — 1. 

The characteristic, or Newton polynomial, of the dimer model on the torus takes the form 

P(z,w)= ^ N„,,„,„ (-l)"^+"-+"^"-z"--w"-, (8) 

where the N„^^„,„ count the number of matchings that have weight (height change) (n^, «„,). 
Furthermore, the partition function which counts the total number of matchings for the 
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Figure 1: Example of a square graph on the torus 

square graph on the torus is given by 

Z = ^{-P{1,1) + P{1,-1) + P{-1,1)+P{-1,-1)), (9) 

where the first term is always zero. P{z,w) evaluated in z, k; = ±1 corresponds to the 
contributions of the four different spin structures. 

3 The dimer model as a free fermion 

As already mentioned in the introduction, many clues point towards the equivalence of the 
dimer model and the free fermion living on the same lattice. On one hand, we know that 
the dimer model corresponds to the two-dimensional Ising model. In the continuum limit, 
the critical Ising model corresponds in turn to the free Majorana fermion. Another clue is 
that the partition function of the dimer model on a genus g Riemann surface is a linear 
combination of Pfaffians, one for each of the different boundary conditions. The one-to- 
one correspondence between the 2^^ equivalence classes of Kasteleyn orientations and the 
spin structures was proved in ||T9l. 

Note that our construction given in the following is different from the transfer matrix 
approach used in Il29ll30l 131)1. 

3.1 The Dirac operator and the Kasteleyn matrix 

In this section, we show that the dimer model on a square lattice can be naturally mapped 
to the dynamics of a free massless fermion on the same lattice. More precisely, we will see 
how introducing a Kasteleyn orientation on the graph can be thought of as the projection of 
the d operator on 7 matrices, see also Il32ll6il. 
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Let us consider the theory for a free massless fermion on a bipartite square graphj^ At 
this stage, the construction is independent of the boundary conditions, i.e. we can think 
of the model as living on the infinite plane. The graph being bipartite, it is natural to use 
staggered fermions [35 1 and associate a two-component real spinor to each pair of nodes: 

Y = (^'^ . (10) 

The naive discretization of the free massless fermion action becomes then 

S= / dxY(x)aY(z)= X];^(x)0(x,y);c(y), (11) 

where the sum runs over all neighbouring vertices x and y on the lattice. The discretized 
Dirac operator can be cast into the form 

a(x,y) = 7'^^^, (12) 
1'^ y\ 

where 7^' are two matrices satisfying the Clifford algebra { 7^^, 7^ } = 2 S>^^ in two Euclidean 
dimensions. On the square lattice, each node has four neighbours, separated by unit vectors, 
so the action reads 

S= E tx{^)^'''\x{x + e'), (13) 

where is the vector e'^ = (cos sin ^). If x is a • site, x ± e'*^ is a o site and vice versa. 
This means that we can split the action into two pieces: 








do. 






Xo 



x»-d»<^-Xo + xo-do.-x»r (14) 



where x» is the vector of Grassmanian variables living on the • nodes. Let us consider the 
first term. The operator 0,o is non-vanishing if and only if • and o are neighbours, and 
takes values proportional to the components of the 7 matrices. More precisely, if we think 
of d,o as a matrix, each line will contain four non-vanishing entries { 7^^, 7^0, —7^0/ — jlo }■ 
Choosing the representation 

these components read { 1, z, — 1, — z }. These are the weights of the four links around each • 



different definition for the Dirac operator, such as the one given in l33ll34l would be needed to deal with 
a general random graph. This goes beyond the scope of the present note. 
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Figure 2: Dirac operator 



node. Since their product equals to (—1) and on this graph the product of weights around a 
node is the same as the product around a plaquette, this is precisely a Kasteleyn orientation 
(see Fig. §3 

The same construction can of course be repeated for the second term in Eq. ( [T4| . In this 
case, we get for each line of the operator do, the four non-vanishing elements {7o,/7o,/ 
— 7o,, — 7o,} which, in the same representation for the 7 matrices that we used above, read 
{ 1, —I, —1, 1 }, yielding again a Kasteleyn orientation. 

Adjacency matrices with equivalent Kasteleyn orientations have the same determinant 
(up to an overall sign), therefore we have 

det{d. o) = det{do.) =Fi{K), (16) 

where X is a Kasteleyn matrix for the graph. Writing the partition function for the free 
massless fermion gives 

Zfermion = J = = ^det(0.o) det(0o.) = Ff{K) , (17) 

which is precisely the number of perfect matchings on the lattice for a given boundary con- 
dition. 

Explicit expressions for both the dimer model on a square graph and the free fermion 

^This argument can be easily generalized to any regular tessellation of the hyperbolic plane (leading to 
higher-genus graphs). One can show that for an n-gon, the product of the weights around a node is (—1) if 
n = mod 4 and (+1) if n = 2 mod 4. For the hexagonal lattice (genus one), the weights around a • site read 
{ i^g27n/3 g-27T!/3 J gj^j their product is (+1). It is worth to remark that for higher genus surfaces, there is no 
unique global choice of complex structure. 
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theory are known. Take a triangular lattice describing a torus with periodicity 



z = z + iN + jMe'^ , 



i,j e z, z e c. 



(18) 



i.e. with modular parameter t = ^e'^. In f36\, it is shown that the partition function for a 
vortex with winding {a, b) on this lattice, which can be mapped to a free fermion, is given by 
the infinite product 



M/2 N/2 

•'^vortex 

*:=-M/2/=-N/2 



a COS 2tz— + 71— + ttcos 2tz— + tz— + 
^ M M/ \ N N ' 



+7 cos 2/1 — - — + TT — - — 



, (19) 



where 



1 — cos ( 

a = i6 = ^ — 

^ sin^ 



_ COS0 

^ ~ sin^e 



2 — cos ( 



sm 



(20) 



Specializing this formula for a square torus with 9 = n/2 such that a = /3 = 1, 7 = 0, = 2, 
we find: 



M/2 N/2 
Zvortex[a(M,N,f) = n n 2 

*:=-M/2/=-N/2 



2kn n\ f 2ln , tt , ^ 



, (21) 



which is precisely the partition function of the dimer model. 



3.2 The scaling limit 

Having shown that the dimer model on a graph is equivalent to a free Majorana fermion the- 
ory we expect the torus partition function to be modular invariant in the thermodynamical 
limit, and to consist of theta functions. 

Since this is a universal property, we specialize to a particular lattice, for example the 
M X N square graph. The partition function for this model was discussed in jJllZl and takes 
the form 



Zmn = \ Yli ^~'^) 

a,b=0 



a+b+ab 



(22) 



where 



Z 



mnU 



M/2 N/2 

n n 2 



2kn 7i\ f2ln , n 

cos I — — + a— + cos — — + 0— 
M M \ N N 



(23) 



J:=-M/2;=-N/2 

If we regard the M x N lattice as a discretization of a square torus with modular parameter 



10 



T = it into M X N squares, the scaling limit is obtained by taking 



M,N ^ CO 



M 
N 



t, t fixed. 



(24) 



This corresponds to sending the lattice size ^ to zero. 

As we have previously pointed out, this is precisely the same expression found for a 
soliton on a triangular lattice [37|. After a careful analysis, the authors prove that in the 
continuum limit, this reproduces the Ray-Singer result for the 9-torsion of the torus 
(see also Il39ll40ll ): 



M,N^co 
M/N=-iT 



t]{T) 



n 2 



(25) 



Instead of reproducing the proof, we will limit ourselves here to a heuristic argument. 
In the infrared limit, the cosine can be approximated by a parabola and the functions above 
behave as follows: 



M/2 N/2 

^mng]= n n 2 

k=-M/2l=-N/2 



2kn n \ f lln , n 

cos I —— + a— + cos — - + 0— 
M My V ^ N 



M/2 N/2 o_2 

n n ^ 

k=-M/2l=-N/2 ^^'^ 



k + 



a 



+ I/ + 2I 



(26) 



In the ^-function regularization, one recognizes the usual product expansion for the theta 
functions. More precisely, introducing t = it we find 



M/N=-n 



1 



2^MN = 2 ^~'^) 
a,b=0 



a+b+ab 



z 



nn 

kGZlGZ 
MNia 



m,n— >oo 2 
m/n=-n a,b=0 



\ E (-1) 



a+b+ab 



?/(t) 



(27) 



ZDirac(T). (28) 



The final result is that in the scaling limit, the partition function is the one of a Dirac fermion, 
or equivalently, for two Majorana fermions. It is not surprising that we encounter a fermion 
doubling effect 14111421 since we have started with a periodic dispersion relation. 



3.3 The loop expansion of the fermion determinant 

After having shown the equivalence of the discretized Dirac operator and the Kasteleyn ori- 
entations and the appearance of the Dirac fermion partition function in the scaling limit, we 
will now directly relate the dimer model to the free fermion via a diagrammatic expansion 
of the fermion determinant. 

Consider a Grassmanian field x living on the nodes V{'i^) of a graph ^ described by a 
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quadratic action: 



(29) 



As we have seen in Eq. (10 we can associate a two-component Majorana fermion Y 



to each black-white pair in a bipartite graph. The path integral for such an action is easily 
computed and is given by the Pfaffian of the matrix A: 



■^fermion 



n d;^. 



£2 



Pf(A), 



(30) 



where x. is the vector x = {XifX2,- ■ }• Iri fact, expanding the exponential (and 

using the fact that Xx are Grassmaruan variables), we find: 



•Zfermion — ^ 



n (1 - XxA^,yXy) 



(31) 



x<y 



All the factors in the product commute, therefore we can order them by increasing x. Ex- 
panding, the terms that survive the integration are those which contain each variable exactly 
once, hence they are products of the | V(^^) | /2 matrix elements Ax,y. There is one such prod- 
uct for each of the permutations of the nodes. Taking into account the signs, we recover the 
expression for the Pfaffian in Eq. (|4]|. 

In the following, we give a diagrammatic expansion for the partition function in terms 
of loops on the graph. Instead of the Pfaffian, we can thanks to relation ([7]| expand the 
determinant under the square root. By definition, the determinant is a sum over all possible 
permutations n of the lattice points, see Eq. (|5]|. The signature sign (tt) of a permutation is 
easily found by expressing n as the composition of cyclic permutations: 



TT = TTs o TTs-i o ■ ■ • o 71/ o ■ • • o K\, 



where tt/ is a cyclic permutation of \r elements. We define 

sign(7r,) = (-1)''+^ 



and 



sign(7r) = nsign(7r/). 

/=i 



(32) 



(33) 



(34) 



Each cyclic permutation n\ can be represented by a closed loop in the graph by joining 
the vertex x by a directed line to the point 7r/(x), see e.g. ||20| . Each vertex in V{^^ can be: 

• isolated if the permutation does not act on it; 

• joined to another one by a double line if the permutation interchanges the two; 
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Figure 3: The diagrammatic representation for the permutation n = o 7x2° ti\ with tti 
(17,18) ^ (18,17), 712 : (8,9,10,13) ^ (9,10,13,8) and 713 : (1,3,11,6) ^ (3,11,6,1). 
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Figure 4: The diagrammatic representation for a permutation allowed by the action for a free 
massless fermion. Orientations are not represented. 



• joined to two others as part of a cyclic chain. 

In particular, a permutation of two neighbouring points is a dimer (with a double line). 

A general permutation n is represented by a graph obtained as the disjoint union of such 
loops, and the determinant is written as the weighted sum of all possible such graphs. Note 
that two points can be joined by a line even if they are not neighbouring on the lattice. If we 
do not impose further constraints on the form of the matrix A, we have a sum over all the 
|y(^^) |! possible permutations and the general term will not be a subgraph of 5f (see Fig. [3]). 

A given graph only contributes to the determinant if the product of the along its 

edges is non-vanishing. This means that imposing constraints on the form of A reduces the 
number of terms needed. 

• A is a local operator if A^^y = when x and y are not neighbouring. In this case, the 
only loops contributing to the determinant are those constructed as unions of edges in 
E{^), and we can limit ourselves to summing over subgraphs of 

• If A does not allow self-interactions, i.e. Ax^x = 0, Vx G V[W), we can consider only 
permutations such that n{x) ^ x, Vx S V{'^). Diagrammatically this means that we 
only consider subgraphs without isolated points. 

The Dirac operator on the graph in Eq. p2| enjoys both properties. In other words, when 
representing the path integral for massless free fermions diagrammatically, we consider all 
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|Y) |m2)) 

Figure 5: Decomposition of a subgraph |Y) containing two loops as the combination | Y) = 
I mi)) (8) |m2)). Orientations are not represented. 

subgraphs containing all the vertices, composed of closed loops or dimers, and without any 
isolated points. An example of such a configuration is depicted in Fig. |4j 

On a bipartite graph, each of these contributing loop configurations can be seen as the 
composition of two perfect matchings. This decomposition is indeed unique. Let |m;)) and 
\mj)) be two perfect matchings. We define the combination |Y) = |m,)) ® |my)) as the col- 
lection of all dimers belonging to either of the matchings, taken with the usual • — > o ori- 
entation if the edge is in |m,)) and with the opposite one if it is in \mj)). Any subgraph |Y) 
contributing to the determinant can always be written in this way. In fact: 

• if |Y) does not contain any loops, its edges are in one-to-one correspondence with 
those of a perfect matching \m)) and we define |Y) = \m) = \m)) (8 |m)); 

• if |Y) contains one loop, this loop consists of an even number of edges since the graph 
is bipartite. Therefore we can decompose it into two disjoint sets of edges that join all 
the vertices. It follows that we can always find two matchings |m,)) and jmy)) that only 
differ by a shift by one dimer on the loop and are equal on all the isolated dimers in | Y) . 
Note that the state | Y) containing the same loop with reversed orientation corresponds 
simply to |Y) = \mj)) (g) |m,)); 

• if |Y) contains n loops (which we regard here as unoriented), it can be constructed as 
the composition of any of the perfect matchings containing the disjoint sets of dimers 
making up each loop. It is easy to convince oneself that there are 2" such combinations, 
corresponding to the 2" possible orientations of the loops (see Fig.|5]). 

Note that the composition of two perfect matchings always generates an allowed config- 
uration. In fact, in each of the matchings, every point in V['^) is touched by an edge and in 
the resulting subgraph, each point is either reached by two distinct lines (and then is part of 
a loop) or by a double one (and is part of a dimer). 

If on a graph ^, there are k perfect matchings, the determinant is written as the sum over 
subgraphs. 
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3.4 Loops as f ermionic states 

As we have seen in the last section, the partition function for the free fermion can be ex- 
panded as a sum over loop configurations on the lattice. Here we show how these configu- 
rations can be understood in terms of fermionic states. 

Consider a square graph on the torus with N horizontal and M vertical cells, each iden- 
tified by a vector z = {z,w) = z ei + w e2- Let h{z) be the height function defined on a 
loop configuration. The function can be labelled by two integers p, q, corresponding to the 
following periodicity conditions: 

/2p,,(z + niNei + naMei) = /2p,(;(z) + nip + n2'?, Mni^X. (35) 

The integers (p, (j) are clearly the winding numbers for the loop configuration!^ 

We would like to consider one of the two coordinates (say z) as Euclidean time and hence 
see the system as the description of a fermion in 1 + 1 Euclidean dimensions. To each pla- 
quette we associate the operators fl(z), fl^(z), h{z), &^(z) living on the left edge, with the 
following interpretations. The operator fl^(z) creates a particle at point w for all times z' > z 
and raises the height function by 1 on the right of point z, while the operator h^{z) creates 
an antiparticle at point z and lowers the height function by 1 on the right. We impose the 
usual anticommutation relations 

{fl(z),fl+(z')} =^z,z', {h{x),h'{z')} = 5,,,,, (37) 

all the other anticommutators vanishing. In particular, we recover the constraint (fl^)^ = 
(b^)^ = 0, which can be translated to the fact that the height function cannot jump by more 
than one unit when passing to an adjacent plaquette. As usual, one can construct a highest- 
weight representation for the algebra by introducing a state |0) which is annihilated by all 
the a and h 

fl(z) |0) = b(z) |0) = 0, (38) 
and defining a general state as the result of the action of a sequence of a' and on |0): 

|Y) = b^(zi)...b^(zNjfl+(zM,+i)...fl'(zM,+Nj|0) . (39) 
^One can define a consistent wavefunction for a Z;^ parafermion as 

^W(,)^,^Mz). (36) 
Here, we concentrate on the k = 2 fermionic case, but the formalism can be extended easily. 
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Figure 6: Height function for a loop state. 



In this operator formalism, the height function is written as 



where we introduced the number operators Nh(z) = fl^(z)fl(z)|^ 

Even though one can reproduce any height function by just using the a and b operators, 
the relations in Eq. ( [37| are not yet enough to describe all the states. This is because we 
have arbitrarily chosen a time direction. One can show that the anticommutation constraints 
should be supplemented by a set of non-local constraints in order to ensure the consistency 
of the height function. Yet it is more convenient to introduce two other sets of operators c, 
and d, which when acting on a point (z, w) create a particle (respectively an antiparticle) 
at fixed time z for each ro' > w. They also satisfy anticommutation relations of the type 



{c(z),c+(z')} = c5. 



{d(z),d+(z')} 



(41) 



and anticommute with all the fls and bs. In this way, we again obtain (c^)^ = (d^)^ = 0, 
which means that the height function cannot jump by more than one unit when passing from 
a plaquette to the one on top or below. Redefining the vacuum |0) as the state annihilated 
by fl(z), b{z), c(z) and d{z), we finally find that a general state is written as 

|Y) = d+(zi) . . . d+(zNjc+(zN,+l) . . . C+(zN,+Njb^(zN^^+N,+l) • • • &^(zn,+N,+nJ 

fl+(zM,+N,+Nb+l) • • • fl^ZN^+N,+N„+Nj |0) • (42) 

Some consistency conditions must be satisfied. In particular, it should be verified that 
the actions of b'^ (z + 1, w) a' (z, w) and of d"'' (z, w + 1) c"*" (z, w) change the height function in 
the same way, namely they raise h{z, w) by one unit (see Figure |6|. It is easy to see that this 

^In a typical Fermi system, h{z, w) is the total charge at point w in space and time z. 
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is equivalent to requiring that on any allowed state |Y): 



h{z,w) |Y) 



|Y) 



Nciz,co)-Ndiz,co) 



|Y) , (43) 



where by allowed we mean that it contributes to the Berezin integral. 

Note that h{z) |0) = 0, and since the operator fl^(z) b^(z) anticommutes with h{z), one can 
easily see that an allowed state |Y) corresponds to a perfect matching if and only if it is 
annihilated by h{z): 

h{z)\PM)=0. (44) 



Other interesting operators can be defined: 

|Y) = [h{N,w)-h{0,zo)] |Y) = 
H,,|Y) = [h{z,M)-h{z,0)]\Y) = 



N 



X]N«(C,a>)-N,,(C,a>) 

M 

J2 Nc{z,a)) -Na{z,w) 



(45a) 
(45b) 



where zv (resp. z) is an arbitrary value. Note that both operators are global since they are 
defined for the whole graph and do not depend on the choice of z or w. We can consider 
h{z) as their local counterpart. Hz, Ha, count the winding numbers of a state |Y) in the z and 
w directions: 



(46a) 
(46b) 



In particular, states with vanishing winding numbers will be annihilated by both Hz and 
Hu,. Another operator is the one that interchanges each fl^ with a b^, each with a and 
viceversa: 



Qi |Y) = Yl \a\z)b{z) + b\z)a{z)] [c+(z)£f(z) + d\z)c{z)] |Y) . 



(47) 



One can easily verify that {Qi,h{z)} = 0. It follows that Qi leaves the perfect matchings 
invariant and pairs all the other states such that Qi \^p,q) = \^p,q) = \^-p,-q)- 

A diagrammatic representation for a general state |Y) is easily obtained as follows. For 
each operator fl"'^(z), one adds a vertical upward arrow on the left edge of the plaquette at 
point z, for each b^(z) a downward arrow on the left edge of the plaquette, for each c"'^(z) a 
left pointing arrow on the bottom of the plaquette and for each d'^{z) a right pointing arrow 
on the bottom of the plaquette (see Fig.|7||. 



In Sec. 3.3 we found the loop states to be combinations of two perfect matchings. We can 
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b^{xi + 1,2:2) 



c\xi,X2) 



Figure 7: Diagrammatic representation for the four creation operators around a plaquette 
with coordinates {z,w). 



indeed recover the same result in terms of this operator formalism, see Appendix |A] for the 
details, where we also give a necessary and sufficient condition for a state |Y) to be allowed. 

Let us summarize the contents of this rather technical section. The terms contributing to 
the expansion of the path integral in Eq. pT) ) can be represented as resulting from the action 
of a set of operators ,h'^ and d"*", defined on the lattice, on a vacuum |0). This set is 
redundant in the sense that the knowledge of the and operators, plus the consistency 
conditions in Eq. (TTTj) uniquely identifies the and . 



We interpret the two directions in the graph as a Euclidean time (z) and a space coordi- 
nate {w). The operator a^{z,w) creates a particle at time z and position w that propagates at 
the same position for all times z < z' < N. In the same way, h^{z,w) creates an antiparticle. 
The height function /z(z, w) computes the local charge (+1 for particles and —1 for antiparti- 
cles) at point w and time z. The charge gradient Hi„ between the points w = Q and if = M at 
fixed time z is conserved during the time evolution. A loop configuration with winding numbers 
{Hz, Hyj) represents the time evolution of a system with initial charge gradient Hi„, whose net charge 
in any fixed point in space w jumps by Hz between z = and z = N. In particular, a (0, H^u) loop 
describes a time evolution where no net charge is created globally, while a perfect match- 
ing corresponds to a system where no net charge is present at the initial time nor is created 
locally. 



4 Categorification of the Newton polynomial 

What follows is inspired by the categorification programme for knot polynomials of Kho- 
vanov ll2Tll22l[231 . The basic idea of his work is that for a knot K, a doubly graded homology 
theory Hij{K) can be constructed, whose graded Euler characteristic with respect to one 
of the gradings yields the knot polynomial in question, e.g. the Jones polynomial, or the 
Alexander polynomial. The homology groups are constructed through a categorification 
process which starts with a state sum representation of the invariant poljmomial. A group 
is constructed for each term in the summation, and differential maps between these groups 
are appropriately defined. Similar categorifications were subsequently performed for other 
invariants with state sums, e.g. various graph polynomials, such as the chromatic polyno- 
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mial [4311441 . and the more general Tutte polynomial Il45l . 

But exactly what is meant by the term categorification? When we categorify, we replace 
set theoretic concepts by category theoretic ones. Sets become categories, functions become 
functors, and equations between functions become natural isomorphisms between functors 
(fulfilling certain relations). The opposite process, decategorification, is much more familiar 
to us. When we decategorify, we are basically throwing away extra information by taking all 
isomorphic objects to be equal. When we say that all vector spaces of the same dimension 
are the same, we decategorify and are left with a set of isomorphism classes. Obviously, the 
reverse process of inventing the missing information is much more difficult. The category 
FinSet of all finite sets for example is a categorification of the set of natural numbers. For an 
easy to read introduction to the concept of categorification, we refer the reader to 1461 . 

After this general digression, we return to our own system of fermion loops. Its partition 
function, the Newton polynomial, is, once more, a state sum. From the point of view of the 
dimer model, its integer coefficients and polynomial nature are obvious. From the point of 
view of the free fermion, this is much less so. Giving this partition function the interpre- 
tation of a graded Euler characteristic explains its integer nature and gives it a geometric 
interpretation, while at the same time suggesting its relation to an index. 

The task at hand is obvious: construct a bigraded co-chain complex from the fermion 
loop states with a differential map preserving the grading, such that the Euler characteristic 
of this co-chain complex yields the Newton polynomial. 

The fact that we can associate four different charges to each loop state living on the torus 
(see Eq. ( |112| ) gives us a certain freedom in the choice of the gradings in the double complex. 
Our emphasis is placed on the matchings, so we choose to sort the loop states into co-chain 
groups according to the maximum winding number of the two constituent matchings in one 
direction (max, Hp|m;))) and to define an internal grading based on the maximum weight in 
the other direction (max, Hq\mi))). This seemingly arbitrary choice just amounts to picking 
a horizontal and a vertical direction in the double complex. The Euler characteristic, for 
example, is not sensitive to this choice. 

When studying the loop states, there is an obvious observation to make: because of the 
two possible orientations of the loops, all states are paired with the exception of the dou- 
ble line perfect matchings. This very much smells of supersymmetric Quantum mechanics 
(SQM), where the supersymmetric ground states are the only ones which need not be paired. 
Indeed, if one takes the total weight of a configuration to be its energy, the perfect matchings, 
having zero winding, become ground states. 

The Witten index, 

dimT^J^o) - dim 7i:Jo) = Tr ( - 1 ) ^, (48) 
counts the number of bosonic minus fermionic ground states and is an invariant of the sys- 
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Figure 8: A loop configuration on the cylinder. In red: closed loops; in blue: double lines. 
This configuration has winding number k — I = +1 



tern. This is, indeed, also what the Newton polynomial is doing! Once we re-interpret our 
fermion loop gas as a SQM system and construct a co-chain complex such that its boundary 
operator can be identified with the Q operator and states with opposite orientations become 
superpartners, the partition function becomes the generalized Witten index of the system. 
This will be done explicitly in the following sections. 



4.1 The dimer model on the cylinder and supersymmetric QM 



In this section, we show how to construct a co-chain complex out of the loop configurations 
of a free fermion on a lattice, and how this system can be mapped to a SQM system in which 
the loop configurations correspond to pure states. For simplicity, we first confine ourselves 
to a square graph embedded on a cylinder where only one non-trivial cycle is present. A 
typical example of such a loop configuration is represented in Figure M 



We have seen in Section 3.3 that a loop configuration can always be decomposed into a 
superposition of two matchings. Let us introduce the following notation: 

• a matching with weight k is written as \k,a)), where a = 1, . . . , Njt labels the matchings 
of the same weight; 

• we label a loop as \k,a;l,b))e with k > I, e = ±1. This corresponds to one of the 
following combinations of matchings: 



\k,a;l,b)^ 


= \k,a))(i 




ifk > /, 


\k,a;l,h)_ 


= \lh))^ 


•) \k,a)) 


if^ > /, 


\k, a-k,h)^ 


= \k,a))(i 


^\k,b)) 


if a > b, 


\k, a;k,h)_ 


= \k,b))(^ 


5 \k,a)) 


if a > b, 


\k, a;k,a)^ 


= \Ka))(^ 


5 \k,a)), 




\k, a; k,a) _ 


= 0. 







(49) 
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Opposite polarizations e correspond to opposite orientations of the same loop config- 
uration. States of the type \k,a;k,a)_^_ (which are in one-to-one correspondence with 
the perfect matchings) are singled out because they are not paired. 

Since each loop configuration with a given orientation can be uniquely decomposed into 
matchings, these quantum numbers can be assigned to a given configuration by inspection. 
The weight of a matching is given by the number of dimers crossing the boundary of the 
fundamental region in positive direction minus the number of dimers crossing the boundary 
of the fundamental region in negative direction (where we choose • — > o as orientation). 
e [k — I) is the total winding number of the loop configuration (which coincides with the 
slope of the height function as defined in Sec. |2|. States of the type \k,a;k,a) are in one-to- 
one correspondence with the perfect matchings. All other states are paired (i.e. there are two 
polarizations). From now on, we will omit to specify the degeneracy quantum numbers and 
write a loop as \k;l)^ (note the semicolon). 

The quantum numbers can be seen as eigenvalues for the following operators: 

H\k;l)^ = {k-l)\k;l)^, (50) 
K\k;l)^ = k\k;l)^, (51) 
U\k-l)^ = e\k-l)^ . (52) 

They commute on the above loop configurations by definition. It follows that H (which we 
want to identify with the energy) commutes also with the operator 

F = 2K-^(n-l), (53) 

which we identify with the fermion number operator. 

In terms of the loops, the above operators have the following interpretation. H counts 
the (unsigned) total winding number of the state. K counts the maximum attainable total 
monodromy of a loop configuration. The state \k;l)g allows for maximally k overall (i.e. 
positively minus negatively oriented) homologically non-trivial windings, of which k — I are 
actually realized. If H is identified with the energy, k gives the maximum possible energy. 
Each loop with non-trivial winding which is added to the configuration raises the energy of 
the state by one unit. Ground states do not contain non-contractible loops. The matchings, 
in particular, contain no loops at all. The eigenvalue of the polarization operator n is the sign 
of the winding number. The fermion number operator f is constructed such, that positive 
polarization states \k;l)^ have even fermion number 2k, while the fermion number of the 
corresponding negative polarization states \k; I) _ is raised by one, i.e. is 2k + 1. The operator 
(—1) gives a Z2 grading according to the fermion parity. States with positive polarization 
have even parity, while states negative polarization have odd parity. 
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To obtain the structure of a SQM, we need to construct two conserved supercharges Q 
and which satisfy the following algebra (see e.g. Chapter 10 of [|47||): 

{QrQ} = {Q+,Q+}=0 (54) 

{Q,Q'} = 2H (55) 
[F,Q] = Q [F,Q'] = -Q'. (56) 

A possible choice is the following: 



Q \k; l)^ = y/2{k- I) \k; I) _ , | Q+ \k; /)+ = 0, 

Q\k;l)_=Q. \Q'\k;l)_ = ^l{k^\kl) 



(57) 



+ • 



It is also useful to define the operator Qi as in Eq. ( |47| l: 

Qi = Q + Q'. (58) 

It reverses the polarization of a state. 



Q^\k-l)^ = ^l{k-l)\k-l)_^, (59) 

and squares to the energy: 

Ql = 2H. (60) 

The space of states can be decomposed into subspaces according to the fermion number. 
In the following, we will concentrate our attention on Q. 

C'' = Span{ \ k,a;l,b)^ \ F \k,a;l,b)^ = p \k,a;l,b)^} , (61a) 

C2'' = Span{ I fc,fl; , (61b) 

C^^+^ = Span { \k,a;lb)_} . (61c) 



We also define 



One verifies easily that 



C+ = 0C2^ (62a) 

k 

C- = C^^^+i , (62b) 

k 

C* = C^ ®C- =^CP . (62c) 

P 



Q : CP ^ CP+^ (63) 
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(and conversely : C — > CP ^). We can now construct the co-chain complex 

^ C° ^ ^ ... ^ C~ ^0. (64) 

Let us now study the cohomology of such a complex. The 0-energy states \k,a;k,b)^ 
and the negative polarization states \k,a;l,b)_ are Q-closed. On the other hand, negative 
polarization states (with k T) are also Q-exact, such that the cohomology coincides with 
the set of supersymmetric ground states: 

0HP(Q) = Span{ | J H |/c, a; = 0} = Span { \ k,a;k,b)J. (65) 

P 

Moreover, all the states with a ^ h are paired (opposite parities). This means that when 
taking the Euler number, only the \k,a;k,a) states (perfect matchings) contribute: 

^ = Tr{-lf = X:(-l)Mim[HP(Q)] = ^Np. (66) 
P P 

As we remarked above, diagrammatically, ground states are configurations with vanishing 
total winding number. The ones contributing to the Euler number are only those that do not 
contain any loops. 

A different invariant can be constructed by taking the trace of the (— l)^z^ operator 
(which is still a conserved quantity since it is a function of conserved operators). It is easy to 
see that it is given by 

;^(z)=Tr[(-l)^z^j =Y,{-lY z'^ din,[HP{Q)] =Y,z^Nk, (67) 

k,e k 



where p = 2k — ^{e — 1). This is the Poincare polynomial for the sequence in Eq. (64 1 and in 
particular, xi'^) = X- Moreover, x{~^) is by construction the partition function for the free 
Majorana fermion on the lattice. 

It is worth stressing that even if, as it is often the case, the zero-modes alone are enough 
to describe important features of the physics of the system such as the partition function, 
we have complete control over all the states in the theory. In Appendix |BJ the generating 
functions capturing all loop states are given. 



4.2 Example: One square on the cylinder 

To illustrate the above, we now present the smallest possible example on the cylinder, con- 
sisting only of one square. The five possible matchings and their quantum numbers are 
shown in Figure |9j They combine into 25 loop configurations, which are, sorted into their 
co-chain groups according to total winding number, depicted in Table IT] 
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—o • 

|1,2> |1,3> 



Q— 



-O 



|0,1> |1,1> |1,2> |1,3> |2,1> 

Figure 9: Matchings for the square on the cylinder 



The partition function can easily be found by inspection and reads: 

P(z)=z(^-^+3-z^ • 
The generating function for this example reads 

G(q, z, y) = 1 + 3yh + /z^ + (1 + y) [si/^z + 3q (yh + /z^) + (?^/z^ 



(68) 



(69) 



4.3 Loops as operators 

In this section, we introduce a different notation for the loops that will allow us to make 
contact with the geometric and gauge theory descriptions in Sec. |5] and |6] Each perfect 
matching can be seen as a pure state in a quantum system obeying a normalization con- 
dition {k,a\l,b) = Ski^ab- The general (quantum) configuration for the system is obtained as 
a linear combination 

W=LL^k,a\Ka) , (70) 

2 

where Aj. „ are complex parameters normalized as usual Y^j^^^ | A)t „ | =1. 

Using this point of view, a loop can be identified with a map going from a state \k,a) to a 
state \ l, h) and as such is represented in bracket notation by \l, h){k,a \ acting as 



{\lh){Ka\)\k\a')=5i,k,5aa'\hh) . 



(71) 



The linear combination of two loops is still well defined, but a new operation is naturally 
defined by the product of two loops: 



{\l,h){k,a\) X [\l',h'){k',a'\) ^ \l,h){k',a'\5u'5ai,'. 



(72) 



Note that this is the usual path algebra on a complete graph whose nodes are labelled by 
\k,a). It is useful to introduce a new representation for our states. Consider a graph with 
N nodes arranged on the lattice at coordinates {k,a), and represent the loop \l,h){k,a\ as 
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Table 1: Co-chain groups for the square on the cylinder 



the arrow from {k,a) to {l,b). Our system of loops corresponds to a complete graph with 
double lines pointing in opposite directions plus a loop on top of each vertex (this is a state 



graph, see Fig. 10 a)). The co-chain groups (Abelian, freely generated by the loops) receive 
the following interpretation: 
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• C^P is the group generated by the loops mapping from a state with winding p to a 
state with winding smaller than p (arrows pointing left), from a state with winding p 
to another one with a smaller degeneracy number (arrows pointing down), or from a 
state to itself (loops). Fig.[lO|b); 

• C^P+^ is the group generated by the loops mapping from a state with winding smaller 
than p to one with winding p (arrows pointing right), or from a state with winding p to 
another with the same winding and higher degeneracy number (arrows pointing up). 
Fig.llOtc). 



In terms of the state graph, this means that the differential map Q2p '■ C^P — > c'^p+'^ maps 
every arrow with a non-vanishing left component to its opposite and annihilates all the 
others: 

Q2p:\l,b){p,a\^ \p,a){l,b\ , (73) 
Q2p:\p,b){P'a\^0. (74) 

On the other hand, Q2p+i : C^P+i C^P+^ annihilates every map 

Q2p+i : |p,fl)(/,&| ^0. (75) 

It is immediate to see that when considering the cohomology one finds: 

• H^P C C^P is generated by downwards pointing arrows and loops. Fig. 10 'd); 



H^P+i c. C'^P^^ is generated by upwards pointing arrows. Fig. 10 e). 



This implies in turn that taking the Euler character, only the loops survive since upward and 
downward pointing arrows are paired and counted with opposite signs (Fig.[lO|^f)). They are 
in one-to-one correspondence with the nodes and all counted with (+) sign, such that 

X = E(-l)''dim(HP) =J^Np = N. (76) 



4.4 The dimer model on the torus and supersymmetric QM 

We are now ready to pass to the actual case of interest, i.e. to extend the categorification 
construction of Sec. 4.1 to graphs embedded on a torus. Now we have two non-trivial cycles, 
so each matching has two weights, kz and kiv, as explained in Sec. |2] These weights will be 
treated on a different footing, but the final result will restore the expected symmetry. 
The notation is modified as follows: 
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Figure 10: State graph for the dimer model in Sec. 4.2 
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• a perfect matching is labelled as |k,fl)), where k is the vector k = {kz,kz„) and a = 

• a loop is the combination of two matchings and is labelled by |k, a; I, b)^, where kz > h, 
e = ±1. 



k,a;l,&)^ 


= |k,fl))$ 


^\lb)) 


ifk ^ I 


k,a;l,b)_ 


= \hb))(g 


\k,a)) 


if k ^ 1, 


k, a; k,b)^ 


= |k,fl))(^ 




if a > b, 


k, a; k,b)_ 


= \Kb))(i 


5 \k,a)) 


\{a > b, 


k, a; k, fl)_|_ 


= |k,fl))$ 


5 |k,fl)). 




k, a; k,a)_ 


= 0, 







(77) 



where k ^ 1 if fc^ > /z or if { fc^ = Iz, k-^o > } . This corresponds to lexicographic order. 
Again, the loops can be seen as eigenstates for the following operators: 



H|k;l), = (fc,-/,)|k;l),, (78) 

X,|k;l), = fc,|k;l), , (79) 

|k;l), = max{huJzu) |k;l), , (80) 

n|k;l), = e|k;l), . (81) 



As opposed to the example on the cylinder, a loop state now comes with two natural signs, 
ez = sign {kz — h), and ew = sign {k-w — Uo)- For the grading by the fermion number, we 
assign only one sign e to each loop state according to Table |2] This is achieved by 



ez 




e 


+ 


+ 


+ 








+ 




+ 




+ 





Table 2: Sign table for torus loop states 

e = + (1 - el) e-,„. (82) 

By choosing the right order on the perfect matchings of the same weight and choosing a sign 
for some loop configurations of winding number 0, the signs of the overall winding numbers 
around the z and w cycles can be brought to correspondence with the signs e^, ew defined 
above. We assign a sign to an oriented loop by choosing a right handed (counterclockwise) 
system to have positive sign and a left handed (clockwise) to have negative sign. The signs 
of a loop configuration are determined by the overall winding number around the cycles. 



28 



Only in cases with an equal number of loops winding a cycle in opposite directions, a sign 
has to be chosen arbitrarily. 

The fermion number operator only depends on Kz and is defined by 

F = 2K,-^(n-l) . (83) 

States can be again collected into subspaces according to the eigenvalues of F, but in 
this case, we have an extra grading given by A;^,. Since [F, = 0, we obtain a direct sum 
structure, and define: 

M 

CP = 0CP'^ (84) 

CP''? = Span{ |k;l)Jf |k;l), = p,K^|k;l), = ^?} . (85) 

The differential operators (or supercharges) are defined precisely as in Eq. ( [57| 

Q|k;l)^ = V2(Az-/z)|k;l)_, iQ'\k;\)^=0, ^^^^ 

_Q|k;l)_ =0. \QMk;l)_ = y^l {h - 1,) \k;l) ^ , 

and it is immediate to verify that they respect the grading since [K-,a, Q] = [K^v, Q^] = 0. 
Again we concentrate on Q. It can be decomposed into the sum 

M 

Qp = EQp,ci' (87) 

where 

Qp : CP'^ ^ CP+''^ , q = 0,...,M. (88) 
As a result, the C* complex is decomposed into C*^q complexes: 

^cO''? ^C^'^ ^^C~w ^0, (89) 



0. 



The analysis of the cohomology is the same as in the case of the cylinder and one can 
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easily convince oneself that it consists of states of zero energy: 



HP'"? = Span { \k,a;l,b)Jh = h}. (90) 

Once more, only the states that are in one-to-one correspondence with the perfect matchings 
[i.e. the ones that do not contain non-trivial closed loops) are not paired and thus survive 
the projection by the (—1)^ operator: 

X{w) = Tr U-lYw^-"'] = Y^wmp,^, (91) 

where we introduced the variable w to keep track of the internal grading. The partition func- 
tion described in Eq. (jsjl is obtained as the trace of the operator ( — 1)^+^2^!" (— z)*^^ (^—■m'^^-^"- 

P{z,w) = Tr [(-1)^+^2^- (^-^f--] = J^{-l)P+^+P^zPw^Np,^. (92) 

P/? 

Note the presence of the factor (— 9^ which coincides with the signature of the spin 
structure. Because of this, the expression differs from the Poincare polynomial for the graded 
complex 

x{z,iv) = Y^zPw'iNp,^. (93) 

The two polynomials, on the other hand, contain the same information and one can pass 
from one to the other using the identity 

2x{z,zv) = -P{z,w) + P{-z,zv) + P{z,-zv) + P{-z,-w) . (94) 

In particular, xi''^) = ^) ^rid 1) is the overall number of perfect matchings on the 

graph. As already remarked at the beginning of this section, it is worth emphasizing that 
the choice of the cycles is arbitrary. This means in particular that one could have used K^v to 
construct the Fermion number operator and Kz as an internal grading without affecting the 
final result. 

The generalization to graphs embedded on higher-genus Riemann surfaces is straight- 
forward and leads to a (2^ — l)-graded complex. 

A simple example to illustrate the above is given in Appendix [C| 



5 Gauge theory interpretation: Preliminaries 

As mentioned in the introduction, there exists a correspondence between quiver gauge the- 
ories describing D-branes probing singular toric surfaces and the dimer model llMllTSll . We 
will make use of it to re-interpret the loops in terms of the quiver gauge theory and the toric 
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geometry. We will quickly summarize the necessary knowledge of quiver gauge theories 
and this correspondence in the following. 

5.1 The quiver gauge theory 

Consider D3-branes in type IIB string theory probing a Calabi-Yau X. They correspond to 
BPS B-type branes given by points on X. If the D-branes are placed at a singularity, they 
are expected to decay into a collection of so-called fractional branes. The resulting world- 
volume gauge theory can be summarized in a quiver graph Qx as follows. Each constituent 
brane appearing with multiplicity N, corresponds to a U{Ni) gauge group which is repre- 
sented as a node in the quiver graph. The massless open strings stretching between the 
fractional branes correspond to chiral fields in the (N,, Nj) representation and are depicted 
as arrows pointing from the L[(N,)-node to the LJ(Ny)-node in the diagram. For the the- 
ory to be anomaly free, the number of arrows fl,i- going from node to node Nj^ (in our 
convention flj., = — must fulfill 

Y^Niaik = 0, \Jk. (95) 

If the ranks of all gauge groups are equal, this reduces to the number of incoming and out- 
going arrows being equal at each node. 

We will be concentrating on the case of the Calabi-Yau X being the cone over a singular 
toric surface S. 

To fully capture the physics of the D-branes, they must be described in terms of the 
(bounded) derived category of coherent sheaves (X). A full explanation of the machinery 
of the derived category is beyond the scope of this paper, the reader is therefore referred 
to ||48| . Luckily, we can avoid working directly in D^(X) by using so-called exceptional col- 
lections of sheaves supported on a partial resolution of the singular surface S. In our case, 
these sheaves can be mostly thought of as line bundles. The exceptional collections form a 
convenient basis for the fractional branes and it is possible to construct the quiver from a 
given exceptional collection and vice versa [49 1. The whole treatment is based on the fact 
that the derived category of coherent sheaves is equivalent to the derived category of quiver 
representations. 

The open strings in the topological sector between two B-branes A and B are parame- 
trized by ©pExt''(A, B). The mass of such an open string between two branes with the same 
grade is given in string units by 

m^ = l{p-l). (96) 

From this, we learn that the massless open strings (i.e. the arrows in the quiver) are counted 
by Ext^, while Ext*^ counts the tachyons. For p > 1, the strings are very massive and are not 
seen by the quiver gauge theory. 
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If oriented loops appear in a quiver (which is always the case for our quiver gauge the- 
ories because of the anomaly-freedom condition in Eq. (|95|), the quiver representations 
become infinite dimensional, causing the methods employed in [49 J to fail. Fortunately, 
this problem can be evaded by considering acyclic subquivers. Some of the arrows of the 
gauge theory quiver Qx are linked to the properties of the surface S itself, while the others 
parametrize the embedding of S in the Calabi-Yau X. Deleting these latter arrows removes 
all oriented loops. The resulting quiver, Qs, allows an ordering relation and is also called the 
Beilinson or Bondal quiver in the literature. 

The derived category of coherent sheaves on S, D^(S), is equivalent to the derived cate- 
gory of quiver representations for Qs, denoted by ( A — mod) (where A is the path algebra 
of Qs). The basic representations L, of the path algebra A are associated to the fractional 
branes at the nodes of the quiver. The arrows in the quiver are thus associated with the 
Ext^ (L;,Ly). The Ext^ (L„Ly) correspond to relations in the quiver which the paths must 
obey. There is another quiver representation given by the projective objects P„ which are 
the subspaces of A generated by all paths starting at node i. The sheaves in the exceptional 
collections we consider are the projective objects Pi, which are in some sense dual to the 
fractional branes L,. 

Let us illuminate a bit the connection between the quivers Qs and Qx and its geometrical 
meaning. Consider the embedding of the surface S into X via i : S ^ X. A D-brane in S can 
obviously also be regarded as a D-brane in X, i.e. we can map objects in D^(S) injectively 
to objects in D"(X) using i^,. Since there might be more morphisms in D^{X) than D\S), 
some of the open strings between two branes in S may live in X but outside of S. This fact is 
captured by the relation 

Ext\{uLi, ULj) = Extlm, Lj) e Ext| (Ly, L,). (97) 

Thus, the decomposition of the arrows of the full quiver Qx into Ext^ and Ext'^ on S accounts 
for the embedding of S in X. 



5.2 Exceptional collections and helices 

In this section, we collect the definitions on exceptional collections necessary for later chap- 
ters. An (ordered) collection of sheaves {JFq, . . . , jF„_i} on S is called exceptional if 

Ext^(.^„^0=C (98a) 
Ext^(J^i,^,) = 0, p>l (98b) 
ExtliJ^i, Tj) = 0, i>j. (98c) 
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The collection is called strongly exceptional, if Extg{!Fi,!Fj) = for p ^ and complete if 
it generates D^(S). Since the D-branes on S can wrap 0—, 2— and 4— cycles, a collection 
must contain n sheaves to be complete, where n =sum over all Betti numbers (= ;c(S)). 
For physics, the collection being strong means that the basis of fractional branes it generates 
does not contain ghost matter. 

Exceptional collections can be transformed into new exceptional collections by left and 
right mutations, which represent an action of the braid group on the set of possible collec- 
tions. On a neighboring pair of sheaves in an exceptional collection, the mutations act as 



The sheaves Le^E/^i and Re^^jE,- are defined via short exact sequences, see IISOllSH . Left and 
right mutations can be seen as braiding and unbraiding. A helix Ti is an infinite collection of 
coherent sheaves {J^;};gz such that 

(a) for any i S Z, . . . , is an exceptional collection. 



The name stems from the fact that after moving n steps to the right, one is back at the orig- 
inal place up to a translation, n is called the period n of H. Each collection of the form 
{J^i+i, . . . , J^i+n} is called a foundation of the helix and Ti is uniquely determined by each of 
its foundations. 

5.3 The correspondence 

Here we quickly summarize the correspondence worked out in ||T4l[T5| . 

1. Start with a toric quiver gauge theory, given by a quiver graph and its tree-level su- 
perpotential. 

2. Draw the quiver diagram on a torus, such that each of the terms of the superpotential 
corresponds to a plaquette. The resulting graph is called the periodic quiver. 

3. Take the graph dual of this graph (i.e. vertices become faces, faces become vertices 
and edges remain edges). Colour the vertices coming from the faces associated to a 
superpotential term with negative sign black; the ones coming from a positive term 
colour white. The resulting graph is bipartite and also lives on a torus. 

4. Solve the dimer model on this graph by taking the Pfaffian of the Kasteleyn matrix. 
This yields the Newton polynomial. The associated Newton polygon is given by taking 



Ri : (E,vE,+i) (E,+i,Re,+i£!)- 



(99) 
(100) 



(b) R"-^J^i = Ti 
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the exponents of the monomials in zw to be coordinates in the (z, w)-plane. The result- 
ing polygon is exactly the toric diagram of the surface the quiver gauge theory comes 
from! 



The singular surface S is specified by the corners of the toric diagram. These corner 
vertices each correspond to Weil divisors and therefore to line bundles. The fully resolved 
surface is given by adding all the points lying in the intersection of the convex hull of the 
polygon and the integer Z^-lattice to the toric diagram. Each such vertex corresponds to 
an exceptional divisor which is blown up to resolve the singularity. The surface S can also 
be resolved partially by blowing up only part of the exceptional divisors. Passing from one 
partial resolution S' corresponding to the exceptional divisor E' being blown up to another 
partial resolution S" with E" being blown up corresponds to a sequence of birational trans- 
formations (blow down E' and blow up E" instead). 

The vertices of the toric diagram can be labeled with the multiplicities of the perfect 
matchings with the corresponding weight. There are certain patterns in the multiplicities 
which are worth pointing out: 

• The vertices at the corners of the toric diagram always have multiplicity one, since they 
correspond to the highest or lowest weight states in either z or w, which are unique. 

• The multiplicities on the edges of the graph follow the rule of Pascal's triangle (i.e. 
1, 2, 1 or 1, 3, 3, 1, etc.). Indeed, the polynomial consisting only of the monomials 
associated to the vertices of an edge has the form (z + w)", or can be brought to this 
form by multiplication with a prefactor z"°w"'° or by redefinition of the variables. 

For the multiplicities of the vertices in the interior of the toric diagram, a pattern is much 
less obvious. 

We illustrate the above by the example of one square on the torus, which corresponds to 
the quiver gauge theory on Fq in Figure [TT] 
The tree-level superpotential is given by 

W = — 2 Xuc^CA^AB^BD + 2 XcA^AB^BD^DCf (101) 

the Newton polynomial is 

P(z,w) = zw ( -- - — + 4 - z - w] . (102) 
\ z w J 

6 Gauge theory interpretation: Results 



In Section 4.3 we have seen that our fermion loop configurations can be interpreted as maps 



from one perfect matching to another. Geometrically, a marked difference exists between the 
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Figure 11: Web of correspondences for one square on the torus/ quiver gauge theory on Fq 



inner vertices of a toric diagram and those on the boundary. The inner vertices correspond to 
exceptional divisors and are compact. The vertices on the boundary must be subdivided into 
the ones at the corners of the diagram and the others. The latter are also exceptional divisors 
and are of the form x C with possible blow-ups, i.e. have a compact and a non-compact 
factor. The corner vertices are non-compact and have the form C x C plus possible blow- 
ups. It is therefore not surprising that it makes a difference whether we only map between 
internal matchings or if boundary matchings appear on one side of the map. When the size 
of the graph is very large, i.e. N, M ^ co, maps between internal matchings indeed are the 
generic case. 

It was shown in ||T7| that for each internal matching PM; G {PM'"^}, an acyclic quiver 
QpM, is obtained by deleting all edges which are part of PM, from the dual quiver Qx- We 
can construct an exceptional collection of sheaves for each such Qpm, - We associate Qpm, to 
the quiver Qg, corresponding to the partial resolution S, of S obtained by blowing up the 
divisor to which PM, belongs. The acyclic quiver Qs represents the system as seen by an 
observer living only on S, while from the point of view of the whole Calabi-Yau X, the low- 
energy theory remains unchanged, regardless of the resolution. The different Qs, associated 
to perfect matchings belonging to the same internal point can be seen as different ways the 
gauge theory is embedded in the surface S, C X. The acyclic quivers associated to different 
resolutions S, and Sj also correspond to different decompositions of the arrows in Qx into 



Extg and Ext^, see Eq. (97 1, but this difference should rather be seen as due to the different 
topologies of the partially resolved surfaces S„ Sj the exceptional collections are supported 
on. 
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On the quiver side, a map between internal matchings corresponds thus to a map /y; : 
QSi ^ Qsjf while on the geometry side, it corresponds to mapping from one exceptional 
collection to another. These maps are always graph isomorphisms of the full quiver Qx, in the 
sense that the adjacency remains unchanged and only the labels of the nodes are permuted. 
This means that the physical gauge quiver Qx remains unchanged by the loop maps, i.e. the 
quiver gauge theory is too coarse to see their effect. At the level of the acyclic quivers Qs,, 
the maps correspond to gauge dualities. Geometrically, the decomposition of the arrows of 
Qx into Extg and Extg changes, which corresponds to a change of embedding. 

In the case of a boundary perfect matching, cycles will remain in the reduced quiver 
QpMm which causes the quiver representation to be infinite dimensional and the construction 
of the exceptional collection to fail. This failure of the known methods to extract useful 
information for these cases prevents us from giving a satisfactory interpretation of maps 
involving boundary matchings. 

There is one special case in which the loop states admit a direct interpretation in the 
quiver gauge theory. It is the case of loops arising from the composition of two adjacent 
boundary matchings. It was noticed in [16j that this results in a so-called zig-zag path. A 
zig-zag path is a path which turns alternatingly maximally left and maximally right at the 
vertices. This maximal turning condition results in also the edges in the dual graph form- 
ing a closed path, which is not the case for a general path. A closed zig-zag path on the 
dimer model graph therefore results in a closed zig-zag path in the dual quiver, which 
corresponds to a gauge invariant trace operator. 

6.1 Matchings and acyclic quivers 

The problem of finding a minimal set of arcs in a directed graph upon the deletion of which 
the graph becomes acyclic is well studied in mathematics and computer science ||52| . Such a 
set is called a mmimal feedback arc set (fas). There exists a precise relation between the prob- 
lem of finding all FAS of a digraph and the problem of identifying all the perfect matchings 
in its dual graph. 

Consider a bipartite plane graph with N nodes and its graph dual The one-cycles in 
are generated by the plaquettes { pj which correspond to the vertices of Removing 
the edge shared by the cycles p, and pj breaks both cycles. A minimal FAS is obtained 
as the collection of N/2 edges { e,j,y^, shared by disjoint pairs of plaquettes p,. and p;^. 
In the dual graph ^, this corresponds to selecting a set of edges joining all the nodes and 
touching all of them only once. In other words, the FAS in ?f ' is a perfect matching in The 
situation is different when §^ is embedded on a Riemann surface of genus g > 0, because in 
addition to the one-cycles generated by the plaquettes, there are 2g equivalence classes of 
cycles of non-trivial holonomy. In this case, the winding cycles in are generated by the 
zig-zag paths. It follows that being a perfect matching in ^ is only a necessary condition for 
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(c) Perfect matching PM3 (weight (1,0)) and quiver QpMj 
Figure 12: Example: Three inner perfect matchings and associated acyclic quivers 



a set of edges to be a FAS in It was shown in [T7| that a perfect matching corresponding 
to an internal point in the toric diagram is always a FAS. Removing a boundary matching, 
on the other hand, always preserves at least one zig-zag path. 

As mentioned before, the gauge quiver Qx is the graph dual of the graph on which 
the dimer model lives. The acyclic quiver Qs„ is obtained from Qx by deleting all the edges 
contained in the nth perfect matching PM„ G {PM'"'}. The deleted arrows are represented 
by dashed lines and correspond to relations in the quiver, i.e. Ext^s. Such an acyclic quiver 
allows an ordering of its nodes such that there are no arrows between the nodes i and ; if 
i < j. If Qs has more than one starting or end point, an ordering ambiguity arises. Indeed, 
Qs is in general a partially ordered set which requires two ordering relations, say left to right 
and bottom to top. 

The maps fm,n^ corresponding to (0, 0)-loops turn out to have a particularly simple struc- 
ture: they always correspond to a cyclic permutation of the nodes in one of the two ordering relations 
(i.e. up to re-ordering in the vertical direction). Maps of non-trivial weight correspond to a 
general permutation of the nodes. In particular, they do not preserve either of the ordering 
relations. In the following, we will attempt to give these maps an interpretation as gauge 
dualities of the acyclic quivers Qs, . 



An example of the above maps is given in Figure 12 This example is big enough to 



have more than one inner point for the toric diagram. Fig. 12 (a) shows a perfect match- 
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ing of weight (0, 0) and the corresponding dual quiver with the deleted edges depicted by 



dashed arrows. The map from this quiver to the one shown in Fig. 12 (b) obviously cor- 
responds to the cyclic permutation 72'^''^^ : {F,E, G, A, C, B, D, H) {A, C, B, D, H,F,E,G). 
Graphically PMi can be recovered from PM2 by first moving G in front of all nodes its ar- 
rows are attached to, then E, and finally F. In the process, dashed arrows which change 
their direction from pointing left to pointing right become solid and solid arrows which 
change from pointing left to right become dashed. Notice the ordering ambiguity in the 
vertical direction. We could have ordered Qpmz slso as e.g. {C,A,B,D,F,H,E,G). Now 



consider the map from Qpmj to Qpmj/ shown in Fig. 12 (c). It corresponds to the permuta- 
tion f^Y'^ '■ {F, E, G, A, C, B, D, H) (C, D, H, E, A, B,F,G) which has an obviously more 
complicated structure and does not preserve either ordering relation. 

6.2 Seiberg dualities and graph isomorphisms 

In the following, we would like to connect the graph isomorphisms of Qx arising in the 
maps fm^n to known gauge dualities of the acyclic quivers. Seiberg dualities Il53l , which are 
IR dualities between gauge theories, are obvious candidates. It was first remarked in (541 
and later elaborated upon in ||49l|55l that certain Seiberg dualities are realized as a sequence 
of exceptional mutations. Here, we are only interested in producing quivers which are sub- 
quivers of Qx- In general, this will not be the case for Seiberg dualities because the rank of 
the gauge groups can change, as in Eq. ( 104| |. 



We want to show here that the effect of performing two consecutive Seiberg dualities on 
the same node is a graph isomorphism on the physical quiver Qx and a cyclic permutation 
on the nodes of the quiver Qg with respect to one of the ordering relations. Let /q be a 
starting node in the quiver, Aout the set of endpoints of the arrows starting in z'o, and A^^ 
the set of starting points for the arrows ending in fo- The effect of a Seiberg duality on z'o is 
the following 1561 : io is moved to the left of the nodes in Ao^t, the directions of the arrows 
coming to and going from /q are reversed, and the arrows between Aout and are changed 
to preserve consistency. More precisely, for A G Aout and B G A^^, 



<.=«B., (103) 



where fly is the number of arrows going from node / to node / (with the convention fly, = — fl/y) 
in the initial quiver Qs, and a[- the corresponding in the quiver Qg obtained after the duality. 
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(a) Qs (b) Q's (c) Q's' 

Figure 13: Double Seiberg duality around the node zq. 



The ranks of the gauge groups are changed as follows: 



N' 



if 7 7^ k 



(104) 



'0 • 



In general, N-^ 7^ N/^ and the transformation is not an isomorphism of Qx- Let us now apply 
the transformation again on the same node /q. The arrows joining /q with the other nodes 
have been reversed, so the role of Am and Aout is interchanged: 



^out ~ ^in r 



^in — ^ 



out • 



(105) 



The node z'o is now moved to the left of A^^j = and the arrows are changed as follows: 



7' 

lot 



Ad;- — fl.-^g — flBjo 



'^BA = ^BA - «^oB«'^^0 



(106) 



(Uab — '^ioA'^Bio) ~ ^ioA^Bio — ^BA ■ 



The ranks now read: 



N'' 



for i ^ io 



E «U^B -K= H "B'o^B - E «'oaNa + 



(107) 



BGA' 



BGAu 



AgAo 



where we used the anomaly-cancellation condition around zq in Qx- 

As expected, applying the Seiberg duality twice just amounts to an isomorphism of the 
gauge quiver Q'^ ~ Qx- On the other hand, the node io has been changed from a starting 
point in Qs to an end point in Qs", and continuous and dashed arrows from/ to z'o have been 



interchanged (see Fig. 13 1. 
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Figure 14: Plaquette move in the square graph 



6.3 Interpretation of the (O,0)-loops 
6.3.1 Plaquette moves 

With the results of the last section, we can now attempt an interpretation of the (0, 0)-loops, 
which form the diagonal component of the cohomology H^'^''^' for the complex in Eq. (|89|. 

Consider a loop configuration which contains the boundary of a basic plaquette and 
otherwise only consists of double-line perfect matchings. We define a plaquette move to be 
a map that switches between the two perfect matchings making up the above plaquette 



loop while leaving the rest of the configuration invariant, see Figure 14 This is obviously a 
(0, 0) map and corresponds to the minimal loop. It can be shown that the set of loops with 
windings {p,q) fixed is connected, with the arcs given by the plaquette moves. It follows 
that all maps of weight (0, 0) are generated by plaquette moves, i.e. we can move between 
all perfect matchings of fixed weight (p, q) by a sequence of plaquette moves. 

The incoming and outgoing arrows of a node io in Qx correspond to the edges of the 
plaquette zq in the dual dimer model graph W. From this dual point of view, the double 
Seiberg duality described in the last section is precisely a plaquette move since the deleted 
edges around io are interchanged while all other edges are left untouched [12]. Since any 
(0, 0)-loop can be decomposed into a product of plaquette moves, it can be interpreted as a 
sequence of double Seiberg dualities on the acyclic quiver. 

It was shown in |57| that a cyclic permutation of the nodes of Qx corresponds to choosing 
a different foundation of the same helix. So whenever we map between perfect matchings 
belonging to the same internal point of the toric diagram, we remain in the same helix. 

Moreover, as we already pointed out, a Beilinson quiver Qs C Qx corresponds to the 
embedding of a (partially resolved) singular surface S in the Calabi-Yau three-fold X. We 
can therefore read a (0, 0)-loop as a map between two possible embeddings i : S "-^ X. 

6.3.2 Relation to crystal melting and black holes 

The dimer configurations of the infinite hexagon graph are in one-to-one correspondence to 
the configurations of a melting crystal corner ||8l. It was shown in | [T2| that they also corre- 
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spond to certain BPS black holes given by D-branes wrapping collapsed cycles in / x 
Zi„ in the large n limit. The black holes are generated by non-geometry changing flop tran- 
sitions, which in the dimer picture are identified with so-called local dimer moves. The local 
dimer moves of fV2\ correspond precisely to the plaquette moves of the hexagon graph. The 
S\ transformations on the fractional branes to which their local dimer moves are related, 
correspond to the double Seiberg dualities described in Section 6.2 Thus, when specializ- 
ing to the hexagon graph, the (0, 0)-loops are the non-geometry changing flop transitions 
in C^/Z„ X Z„ which generate those BPS black hole charge configurations that are parame- 
trized by the 3d crystal melting configurations. On a plane graph, everything is captured by 
the (0, 0)-loops since no non-trivial winding states exist. 



6.4 Interpretation of the loops with non-trivial winding 

As mentioned before, the maps corresponding to loops of non-trivial winding result in a 
general permutation of the nodes of the gauge quiver Qx- Geometrically, we map from 
one partial resolution of S to another, which is a birational map. The associated exceptional 
collections produce equivalent derived categories. 

Unlike the (0, 0)-maps, loops with non-trivial winding can in general not be written as a 
sequence of Seiberg dualities. We have seen that plaquette moves cannot generate winding 
loops, therefore we cannot write the latter as a sequence of double Seiberg dualities. Single 
Seiberg dualities in general change the ranks of the gauge groups, see Eq. ( |104[ |, and thus do 
not lead to a graph isomorphism of the gauge quiver Qx- 

Since the full quiver gauge theory remains unchanged, the action of the weight changing 
maps filff' on the acyclic quivers can still be seen as a generalized gauge duality on the 
theory on the (partially resolved) surface S. 

What is the geometrical meaning of this? Morally speaking, an exceptional collection 
forms a basis for the physical branes. All partial resolutions lead to the same quiver gauge 
theory in X, so at the level of the derived category, the associated exceptional collections are 
equivalent. Passing from one to the other corresponds to a change of basis. We have seen 
that with the (0, 0)-maps, we change between foundations of a helix. This is a very particular 
change of basis, which preserves the order of the sheaves in the collection and only amounts 
to tensoring the canonical bundle a certain number of times to all sheaves in the collection. 
We have seen on the other hand, that the general permutation of nodes which takes place 
in a weight-changing map does not preserve the order of the sheaves in the collection. This 
means that we no longer remain in the same helix. At the level of the sheaves, this map is 
realized by a general sequence of mutations. 
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7 Conclusions 



In this note, we relate a statistical mechanical system, the dimer model on a torus, to a sys- 
tem of lattice quantum field theory, i.e. the massless free fermion. The loop configurations 



obtained after a diagrammatic expansion of the fermion determinant (Sec. 3.3 1 can be inter- 
preted as states in SQM after employing categorification techniques inspired by Khovanov's 
work (Sec. |4j|. The Newton polynomial of the dimer model becomes the generalized Eu- 



ler characteristic of a co-chain complex (Sec. 4.4 1. The states with vanishing overall winding 



become the supersymmetric ground states of the system. Since all loop states except the dou- 
ble line perfect matchings are paired, taking the Witten index of this system again counts the 
number of perfect matchings (Eq. ([92|). 

In the following, we make use of the dimer model-quiver gauge theory correspondence 
to give the loop configurations an interpretation in yet another picture. Interpreting the 
loops as maps from one perfect matching to another (Sec. |4.3| , maps between internal perfect 
matchings associated to the same point in the toric diagram become changes of foundation of 
the associated helix and sequences of double Seiberg dualities in the acyclic quiver (Sec. |6.2[ |. 
When specializing to the honeycomb graph, the plaquette moves which generate all (0, 0)- 
maps correspond directly to the BPS black hole configurations which are parametrized by 
crystal melting configurations (Sec. 6.3.2| . 



Mapping between matchings associated to different internal points is a more severe 
change and cannot be achieved as a sequence of Seiberg dualities. In the language of ex- 
ceptional collections, we map from one helix to another. 

In Table |3| we summarize the results of this paper by giving a dictionary between the 
loop states, the fermion states, the gauge theory and the geometry. 

From this point, one can go to several different directions for further research. 

• The dimer model admits a description as a free fermion model based on a transfer 
matrix approach (see l|29ll30ll3T]| and also lU) which, at first sight, is different from the 
one we introduced in this note. It would be very interesting to clarify the link between 
these two descriptions. 

• Furthermore, it would be interesting to use a more general definition for the Dirac 
operator on a graph than the one in Eq. p2| , to prove the equivalence to the dimer 
model completely independently of the type of graph and its representation. Indeed, 
this could provide an elegant solution to the problem of choosing a Dirac operator on 
a random graph that depends only on the the adjacency. 

• The algebraic structure we introduce in Sec. |4] also deserves further investigation. A 
better understanding of the complex and the mapping to SQM might indeed cast a new 
light on the dual quiver gauge theory. 
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• A final word must be spent on the geometric interpretation in Sec. |5j To our knowl- 
edge, the meaning of the boundary perfect matchings (and, correspondingly, of the 
infinite-dimensional representations of the path algebras) in terms of derived category 
is not yet clear. This is an obvious prerequisite to a complete understanding of the 
geometrical meaning of all the loop states that appear in the theory. 

The study of the dimer model, which due to its combinatorial nature provides a relatively 
easy playground, has proven to be very fruitful thanks to its manifold connections to other 
physical systems. We believe that these inter-relations have not yet been fully exploited and 
deserve further attention. 
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A Perfect matchings in the operator formalism 



In Section 3.3 we found the loop states to be combinations of two perfect matchings. It 



would be interesting to recover the same result in terms of the operator formalism of Sec- 



tion 3.4 To do so, we start by distinguishing between even and odd plaquettes. We will say 
that a plaquette with coordinates z = {z,w) is even (resp. odd) if z + if is even (resp. odd). 
We use the definition of the height function from the superposition of two perfect matchings 
PMi — PM2 given in Section |2] If an odd plaquette contains either an upward arrow on the 
left side, a downward arrow on the right side, a right pointing arrow on the bottom side or 
a left pointing arrow on the upper side, it belongs to the first perfect matching PMi, in the 
converse case, it belongs to PM2. For an even plaquette, the arrows are reversed. The strat- 
egy to construct an operator which only preserves the operators corresponding to one of the 
perfect matchings, say PMi, when acting on a state |Y) is simple. Consider e.g. all the black 
nodes. For each of them add all the annihilators corresponding to the arrows belonging to 
PM2. If |Y) is an allowed state, only one operator is paired for each node which gives a 1 as a 
result of the anticommutation, while all the others will just annihilate the state. This results 
in only the operators belonging to PMi surviving and so yielding the result we are looking 
for. We can define the operators PMi and PM2 explicitly as follows: 

N M 



IV Mil j _-\Y^ 

= n n 2 — [^^^' ^ ^ ^^^' ^ - 1) + - 1' ^)] ' (108a) 

Z = l IV = 1 
Z f 1 _|_ ( 1 + 

PM2 = n n 9 \b{z,w)^c{z,w)^a(z,w-\)^d{z-\,w)\ . (108b) 



Z = l IU=1 
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A state |Y) = \nii)) \m2)) is decomposed into its two constituent perfect matchings by their 
action: 

PM, |Y) = |m,)) , i = l,l. (109) 

Now we are able to give a consistency condition for a state to correspond to a loop appearing 
in the expansion of the fermionic action. The idea is that the state should be decomposable 
into two matchings, which both touch each node once. Let us therefore define 



Pi{z,w) 

P2{Z,W) 



^-^^^^[lSla{z,w) + Nd{z,w) + Ni,{z,w-l)+Nc{z-l,w)-l] , (110a) 



[Ni,{z,w) + Nc{z,w) + Na{z,w - 1) + Na{z - l,w) - 1] 



(110b) 



An allowed state |Y) must satisfy the local conditions 

Pi(z)|Y) =P2(z)|Y) =0, Vz. 



(Ill) 



The action of the height operators in Eq. ( [45| is still defined on the perfect matching states 
|Y)). To be consistent with the constructions in Section [2| we introduce a sign to distinguish 
between Imi)) and Iwii)): 



H,\m,)) = {-iy+' 
H.„|m,-)) = (-1)'+! 



N 



X]N„(^,a;)-Nf,(C,a>) 



M 



£ Nc(z,a;) -Nrf(z,a;) 



|m;)) 



(112a) 
(112b) 



We find that each loop can be identified by four charges, corresponding to the eigenvalues 
of the winding operators acting on the two perfect matchings. The relation between the 
winding of a loop and the weights of the constituent perfect matchings is 



H,|Y) = H,|mi))-H,|m2)), 
H,„|Y) = H,„|mi))-H,„|m2)). 



(113a) 
(113b) 



Extending the definition of the height function to perfect matchings requires some care, since 
the equality in Eq. ( [43| does not hold for perfect matchings. A possible way out consists in 
symmetrizing the sum of the two expressions and define 



}i{z,w)\mi)) 



(-1) 



X]N„(^,i(;)-Nb(^,a;)+ ^ Nc(z,a;) - Nrf(z,a;) 

J=0 a;=0 



(114) 



In this way, one recovers once more 



/2(z) |Y) = /z(z)|mi)) - h{z)\m2) 



(115) 
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B Generating functions 



In this section, we give the generating function which contains the full information of the 
fermion loop gas, i.e. all the loop states. The dimer model partition function and the Witten 
index can be recovered from it. 



B.l Generating function for the fermion loops on the cylinder 

We compute the generating functions for the loops defined as 

G{q,z,y) = Tr 



(116) 



which keeps track of all the information. Expanding the trace one finds an explicit expression 
depending on the matching degeneracies Nj-: 



N 



1 N N N 



G{q, z, y) = y^'z^ + (i + y) <| i ^ (N, - 1) y^^z^ + ^ ^' E y^'z^ 



k=0 



k=0 



1=1 k=l 



(117) 

All the physical information about the system is contained in the generating function G. 
For example: 

• the Euler number for the complex in Eq. (|64]| is = Tr (—1)^ = G(l, 1, —1) ; 

• the Poincare polynomial for the complex is x{^) = Tr (—1)^ z'^ = G(l,z, —1) ; 



• the total number of perfect matchings is N = Tr (— 1) 



G(l,-1,-1); 



the multiplicity of weight-A: matchings is N} 



1 3* 



ij=l,z=0,y= — l 

the dimension of the k-th chain group is dimC^ = ^ ^G{q,z,y) 



^ k\ az'f y 
1 d'' 



q=l,z=l,y=0 



the number of loops with winding number his L 



1 d'' 



h — M dql- 



G{q,z,y) 



q=0,z=l,y=l 



B.2 Generating function for the fermion loops on the torus 



Let us now extend the definition of generating function for the loops to graphs embedded 
on the torus. Following the formula for the cylinder in Eq. \[17) we define: 



G(q,z,y) = Tr 



q z y 



(118) 
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where the notation is a shortcut for = q^^cjw"" and = z^~w^^". Expanding in the loop 
basis we find: 



G(q,z,y)= £ N^y''^z^ + ^^l t 

k={0,0) lk=(0,0) 



(Nk-l)y2'^^z''+ 



N 



N 



+ E ^kNu-i + N,^fr„-/„.)%z-yitJ y''^^" ' (119) 



1=(0,0) k=l 

where the primed sum means that the value 1 = (0, 0) is to be omitted. All the physics 

1 

of the system is encoded into G(q, z, y). In particular: 

• the Euler number for the complex in Eq. ( [89| is xi'^) = Tr [( — 
= G{{l,l),{l,w),-l); 

• the Poincare characteristic for the complex is = Tr [(— l)^z''] 
= G((l,l),z,-l); 

• the number of perfect matchings with height k is Nk = ^ 



where we introduced the notation k! = kzlkj^^. , |k| = kz + kj^ , and 



q=(l,l),z=(0,0),y=-l 



• the dimension of the chain group is dim = 

• the number of loops with winding h is = jj, |^ G 

C Example: One square on the torus 



dykzdw^,^ q=(l,l),z=(0,l),y=l 
q=(0,0),z=(l,l),y=l ■ 



To illustrate the construction given in Section 4.4 we now present the smallest possible ex- 



ample on the torus, consisting only of one square. The eight possible matchings and their 
quantum numbers are shown in Figure |9j They combine into 64 loop configurations, which 
are, sorted into their co-chain groups according to total winding number, depicted in Table|4] 
The partition function of the dimer model is the one given in ( 102| , the generating func- 
tion for this example reads 



G(q,z,y) 



w + y^z + Awy^z + w^y^z + wy'^z^ + 



+ (1 + y) [(6wi/^z + 4(?z + ^^u + qzqw) wy'^z + (4(?z„ + ql, + qzqw) w'^y'^z+ 



+ (4(j2 + ql + qzqw) wy^z^ + qzqww'^y^z^ 



(120) 
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(1,0) 



(0,1) 



|(1,0),1> 

(1,1) 



(2,1) 



— r ^? — — — r — r ^ — t 

|(0,1),1> |(1,1),1> |(1,1),2> |(1,1),3> |(1,1),4> |(2,1),1> 

(1,2) 



|(1,2),1> 



Figure 15: Matchings for the square on the torus 
Table 4: Co-chain groups for the square on the torus 





CO 


ci 


C2 


C3 


C4 






|1,2,1;1,0,1), |1,2,1;1,0,1)_ 



|1,2,1;1,1,4), |1,2,1;1,1,4)_ 



2 




ll 2 1-1 


1,3)+ |1,2,1;1, 


1,3)_ 






2 




|1,2,1;1,1,2)+ |1,2,1;1,1,2)_ 







Continued on next page 
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Table 4 - continued from previous page 



K,, C° O 



2 


1,2,1; 1,1,1) + 

n 


1,2,1;1,1,1)_ 

n 


2 


|1,2,1; 1,2,1) 


|2,1,1;1,2,1)+ |2,1,1;1,2,1)_ 


1 


1,1,4;0,1,1) + 


1,1,4;0,1,1)_ 


1 




i;0, 

1 : 


1.1)+ 


1,1,3;0,1,1)_ 


1 


1,1, 2; 0,1,1) + 


1,1,2 


i;0, 

1 N 

)^ M 


i,i)_ 
r 


1 


1,1,] 


;0, 

1 c 

) M 


1.1)+ 
r 


1,1,] 


;0, 

1 ( 

— 1 


i,i)_ 
r 


1 


1,0,1; 0,1,1) + 


1,0,1;0,1,1)_ 


1 


114-1 

-■^z ^/ ^/ -"^z 


0,1) + 


114-1 

-•^z -"^z ^/ -"^z 

H 


o,i)_ 


1 


— 


i; 1,0,1)+ 


1,1,2 

— 


i;i,o,i)_ 


1 


1,1,2; 1,0,1) + 

1i 


1,1,2;1,0,1)_ 

u ■ 



Continued on next page 



50 



Table 4 - continued from previous page 





CO 




C3 


C4 




1 




|1,1,1;1,0,1) + 


|1,1,1;1,0,1)_ 







|1,1,4;1,1,3), |1,1,4;1,1,3)_ 



1 


- 


|1,1,4;1,1,2) + 

- n 


1 1,1, 4; 1,1, 2) 


1 


- 


|1, 1,4; 1,1,1) + 

- ±± 

t i 


|1, 1,4; 1,1,1) 

±± 

* t 


1 


- 


|1,1,3;1,1,2) + 


1 1,1, 3; 1,1, 2) 

tt 


1 


- 


|1,1,3;1,1,1) + 

- n 


|1, 1,3; 1,1,1) 

u 


1 




|1,1,2;1,1,1) + 

- u 


|1, 1,2; 1,1,1) 


1 




1 1,1, 4; 1,1, 4) 




1 




1 1,1, 3; 1,1, 3) 

- u 




1 




ll 1 

_ — ( 

— c 


2;1 

H 
H 


,1,2) 

3— 
» 





|2,1,1;0,1,1)+ |2,1,1;0,1,1)_ 



^^^^ 



|2,1,1;1,0,1)+ |2,1,1;1,0,1)_ 






|2,1,1;1,1,3)+ |2,1,1;1,1,3)_ 

n n 

|2,1,1;1,1,2)+ |2,1,1;1,1,2)_ 




|2,1,1;1,1,1)+ |2,1,1;1,1,1)_ 
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Table 4 - continued from previous page 



|0, 1,1; 0,1,1) 

1 t*i 



|i, 1,1; 1,1,1) 



**? — r- 



|2, 1,1; 2, 1,1) 




1,0,1; 1,0,1) 
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